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ABSTRACT 

Recent observations of heavy elements produced by rapid neutron capture 
(r-process) in the halo have shown a striking and unexpected behavior: within a 
single star, the relative abundances of r-process elements heavier than Eu are the 
same as the same as those of solar system matter, while across stars with similar 
metallicity Fe/H, the r/Fe ratio varies over two orders of magnitude. In this 
paper we present a simple analytic model which describes a star's abundances in 
terms of its "ancestry," i.e., the number of nucleosynthesis events (e.g., supernova 
explosions) which contributed to the star's composition. This model leads to a 
very simple analytic expression for the abundance scatter versus Fe/H, which is 
in good agreement with the data and with more sophisticated numerical models. 
We investigate two classes of scenarios for r-process nucleosynthesis, one in which 
r-process synthesis events occur in only ~ 4% of supernovae but iron synthesis 
is ubiquitous, and one in which iron nucleosynthesis occurs in only about 9% 
of supernovae. (the Wasserburg- Qian model). We find that the predictions 
in these scenarios are similar for [Fe/H] > —2.5, but that these models can 
be readily distinguished observationally by measuring the dispersion in r/Fe at 
[Fe/H] < -3. 

Subject headings: nuclear reactions, nucleosynthesis, abundances 
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1. 



Introduction 



Neutron capture processes dominate the nucleosynthesis of elements beyond the iron 
peak. The physics of two neutron capture mechanisms has long been understood, and the 
astrophysical site for slow neutron capture nucleosynthesis (the 5-process) has been shown to 
be low-mass stars (e.g., Busso, Gallino, & Wasserburg (1999)). However, the site for rapid 
neutron capture nucleosynthesis (the r-process) has not yet been unambiguously identified, 
although it is very hkely connected to massive stars. The production site itself might be found 
in some or all Type II supernova (SN) events (e.g., Woosley, Wilson, Mathews, Hoffman, & 
Meyer (1994)), or in binary neutron star mergers (e.g., Eichler, Livio, Piran, & Schramm 
(1989); Rosswog, Davies, Thielemann, & Piran (2000)). Truran (1981) showed that the halo 
stars are a particularly useful laboratory for study of the r-process, as the observed neutron 
capture abundance patterns in these stars indicate that the s-process component drops out, 
and the r-process dominates, as one looks at stars with [Fe/H] < —2.2 (Burris et al. (2000)). 
Recent years have shown that halo stars indeed give unique insight into the r-process. With 
the advent of high dispersion, high S/N measurements, the abundance observations within 
and among halo stars has led to surprising new discoveries with important consequences for 
theories of the r-process as well as halo formation. 

The halo star with perhaps the most striking r-process composition is CS 22982-052. 
With [Fe/H] = —3.1, this is an ultra-metal-poor star;^ its abundance ratios through the iron 
peak are typical for a halo star However, this star also has been observed in 20 r-process 
elements from barium upwards (the "heavy" r-process), with abundances (Sneden et al. 



The element-to-element scatter in this ratio is consistent with the observational errors. This 
remarkable trend does not appear to hold for the "light" r-process, 40 < Z < 56 (Sneden et al. 
(2000)). Nevertheless, the fact that the heavy r-process elements agree so well with each other 
is stunning, given that this star could well have abundances that reflect the nucleosynthesis of 
a single supernova, while the solar abundances average over many generations of supernovae. 
The strong implication is that there is a unique astrophysical site for the r-process (at least 
for the heavy r-elements, which for brevity we will refer to simply as "r-elements" hereafter). 

This result has profound implications. First, the mere presence of r-process in very 
metal-poor stars points to origins associated with massive stars. Moreover, given the appar- 



^In the standard notation, [A/i?] = log(^/-B)* — \og{A/ B)q; we take ultra-metal-poor stars to be those 
with [Fe/H] < -2.5. 



(1996; 2000)) 
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ent universality of the r-process, one might expect that all ratios among r-process nuclides 
to always be constant (with pollution by the s-process increasing with [Fe/H]). This leads to 
the prediction that at the lowest metallicities, when the s-process has not turned on, ratios 
among r-process elements should be fixed and the same across all stars. So far this appears to 
be true in the other ultra-metal-poor star for which there is similarly good data, HD 115444 
(Westin, Sneden, Gustafsson, & Cowan (2000)), which has [Fe/H] = -3.0 and [r/Fe] = 0.96. 
Very recently, Cayrel et al. (2001) have shown that CS 31082-001 (with [Fe/H] = -2.91) 
has Os/Fe and Ir/Fe abundances which give [r/Fe] = 1.98 and it! = 96. This would be 
the highest r-process overabundance yet observed. Clearly, a systematic study of r-process 
abundance patterns in this star is of the highest priority (particularly since this is the first 
star in which uranium has been detected). 

CS 22892-052 and HD 115444 have iron abundances that are identical within errors, and 
both show a remarkable constancy in their r/Fe ratios for different r-elements within each 
star. However, these stars show a significant difference in their mean r/Fe ratios. This is an 
example of the large scatter in r/Fe ratios that has been observed in halo stars. As we will see, 
the dispersion in r/Fe ratios is small at [Fe/H] > — 1, but grows with decreasing metallicity, 
finally spanning two orders of magnitude for ultra-metal-poor stars. The observed large 
scatter in r/Fe in Pop 11 stars — and the very high r/Fe in a few — demands that not all 
stars make both the r-process and iron. 

Models for r-process nucleosynthesis and for Pop 11 chemical evolution must account at 
once for the remarkable constancy of r/Fc among different elements within individual stars, 
as well as for the variability of r/Fe between stars. These facts together suggest the following 
simple picture for Pop II r/Fe. We assume that the observed Pop II r/Fe abundances in 
each star reflect r-process contributions of a few supernovae. The basic idea is that 

• Different r /Fe arise from different mixing between the dominant r-process sources and 
the dominant Fe sources; and 

• The r/Fe scatter then reflects the amount of mixing between the two sources. 

The key is the inhomogeneity of the halo, which this scheme in fact quantifies. 
This scenario makes quantitative predictions: 

1. There is a maximum to r/Fc, which is the r/Fc production ratio in those stars which 
produce r-process elements, undiluted by any (r-poor) Fe events 

2. The minimum to r/Fe is the r/Fe ratio in the r-poor events. This minimum may be 
as small as r/Fe = 0. 
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3. The placement of a star between these extremes quantifies a admixture of nucleosyn- 
thesis sources recorded in the star. Thus, (a) at the very earliest times, when single 
SN events are really all that contaminate a given star, the r/Fe scatter should only 
populate the extremes, (b) At later times, but before mixing is efficient, the regions 
between the extremes will be filled in. (c) At still later times, mixing becomes efficient, 
and then the scatter decreases towards a " universal" mix of the two sources. 

If this fine of reasoning is correct, it has the very significant implication that r/Fe becomes 
a tracer of the inhomogeneity of the halo. One can see the transition from single to multiple 
events, and also the transition from multiple sources to a well mixed Galaxy. The details of 
these transitions allow one to correlate metallicity and mixing. And we can at least grossly 
confirm the time sequence of the nucleosynthesis events. 

Several groups have used similar arguments to motivate detailed models which explain 
r/Fe scatter in terms of r-process nucleosynthesis and an inhomogeneous chemical evolution 
of the Galactic halo. Ishimaru & Wanajo (1999) construct a one-zone model for the halo, 
but introduce abundance dispersion through a Monte Carlo realization of individual super- 
nova events with yields of Fe and Eu which are strong function of mass. They obtain good 
fits to the observed scatter for models in which the r-process is produced by a small frac- 
tion of supernovae. Argast, Samland, Gerhard, & Thielemann (2000) use similar r-process 
progenitor masses, but drop the assumption of one-zone evolution. They model the spatial 
inhomogeneity and incomplete mixing of expanding supernova remnants in the halo, and 
again find good agreement with the observed Eu/Fe dispersion. 

The goal of the present paper is to aid in the understanding of these detailed models, 
and to help focus attention on the key physics relevant for this problem. We will show (§2) 
that the r/Fe scatter follows quite generally from the stochastic behavior of the halo chemical 
evolution, and that one can derive simple, explicit, analytical expressions for the abundance 
scatter as a function of metallicity. We will do this for two scenarios, one in which the 
r-process occurs in only a small fraction of supernovae but Fe is ubiquitous (§2.1), and the 
scenario of Wasserburg & Qian (2000), in which the converse holds (§2.2). Fortunately, these 
scenarios can be distinguished by future measurements of r/Fc scatter in ultra- metal-poor 
stars, and we present observational strategies to do this (§3). Discussion and conclusions 
appear in §4. 
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2. Two Scenarios for r^-process Origin 



The r/Fe scatter demands that the r-process and iron are not co-produced in the same 
ratios in the same events, but rather there will be r-rich and r-poor events. Also, very high 
and very low r/Fe ratios are observed. The highest r/Fe event observed places a lower limit 
on the r/Fe yields of the r-rich events. Furthermore, the fact that the highest r/Fe is much 
larger than the Pop II mean implies that the r-rich, high r/Fe events occur at a very different 
rate than low r/Fe events. One can envision either "iron-dominated" schemes in which high 
r/Fe events are rare, but iron-producing events are not; or, alternatively, the possibility that 
r-process events are the more common, with iron-producing events rarer. We now consider 
each of these scenarios in turn. 



We assume (1) that the supernova r-process yields are bimodal, i.e., that some super- 
novae have a high r-process yield and others a low one (possibly zero); and (2) for simplicity, 
we assume not only that all supernova make iron, but also that the Fe yield takes the same 
constant value for all supernovae.^ With these assumptions, supernova ejecta can take one 
of two values of r/Fe. We will denote the objects with the larger r/Fe ratio as class A, and 
those with the lower ratio as class B. We purposely avoid a more descriptive nomenclature- 
e.g., "high" and "low" -as we wish to avoid both confusion (high yields? rates? masses?) and 
prejudice with regard to the nature of the progenitors which are the dominant r producers. 
To simplify notation, we will define the scaled r/Fe ratio to be 



which implies that [r/Fe] = logR. 

Thus, a supernova will produce a fixed amount of iron, and yields r-elements at either 
the high value Ra or the low value Rb- Denote the fraction of class A supernovae (averaged 
over the halo IMF) to be /„; then the fraction of class B supernovae is of course fb — ^ — fa- 

The r/Fe data show that at a metaUicity Fe/H ~ —1, [r/Fe] converges to a value of 
~ 0.3, i.e., R ~ 2.^ We take this value as the mean R attained when averaging over a large 



^It is possible that some supernova make neither r nor Fe; if such events exist, they are irrelevant to the 
present discussion. 

^Above this metaUicity r/Fe decreases, but this is presumably due to the addition of the Fe yields of 



2.1. 



Iron-dominated Models 



R = 



r/Fe 



(2) 
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sample of progenitors, whose fractions reflect the underlying and fb dictated by the IMF. 
Thus we have that 

R — fa Ra + fb Rb — Rb + fa {Ra — Rb) (3) 

We can use the observed scatter to estimate the values of Ra and Rb- The highest observed is 
[r/Fe] ^ 1.7 (eq. 1); we take this as a measure of Ra — 50. On the other hand, the lower limit 
of the r /Fe scatter is not as well determined. The available data suggests that the minimum, 
if any, differs depending on the elemental r-tracer used. At most, we have [r/FcJmin ~ —0.6, 
giving Rfj ~ 0.25; it is also possible that [r/FcJmin ~ —1-7, giving Ri, = 0.02 ~ 0. Upper 
limits to r/Fe have not been reported, but this may be a selection effect. A key issue is 
whether there could have been lower r/Fe values than the ones reported; such information 
places useful constraints on r-process production. For our purposes, as long as Rb <^ R, then 
we may take Rb — 0, as we will now see. 

Once one establishes the observed values for R, Ra, and Rb, then the relative numbers 
of class A and B supernovae is now fixed. Eq. (3) gives 

f _ R ~ Rb ^ R (A\ 

^a ~ -Kb n-a 

where the last expression uses the observed fact that Rb/R <^ 1. For R — 2 and Ra — 50, 

cq. (4) gives fa = 3.5% for Rb = 0.25, and fa = 4.0% for Rb = 0. Thus we see that observed 
fact of the very large r/Fe variations immediately implies that class A events are required 
to be uncommon, regardless of their physical nature. That is, the events which produce the 
Y-process originate from a small fraction, ~ 4%, of all massive stars. 



2.1.1. Basic Model 

Using these numbers, we are now in the position to model the scatter of r/Fe. To do 
this, we create a set of halo stars and deduce the history-the nucleosynthetic ancestry-of 
each. We assume that each of the halo stars we create incorporates gas which has been 
enriched by some number A^"sn of supernova of either type: this is the number of supernova 
ancestors for the star. We create stars by allow iVsN to run from 1 to N-^^, and sample 
equal intervals in A^sn- Note that we are free to chose the number of halo stars we create, 
and that this number is unrelated to the number of supernova ancestors assigned to a given 
halo star. 



Type la supernovae. We do not include these events (and implicitly assume their r-process yield to be zero); 
thus we will not attempt to model r/Fe at [Fe/H] > — 1. 
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The metallicity follows from the assumption that all stars have the same iron yield, so 
that Fe/H oc A^^sn- We introduce a parameter (Fe/H)inin, the iron enrichment from a single 
star. Then Fe/H = Ns^{Fe/RU^, and [Fe/H] = [Fe/RU^ + logTVs^. 

The r-process pattern of the stars is also fixed in terms of the stellar ancestry. We assume 
that the mean value of it! for our population of stars is always We then introduce scatter 
around R by sampling stars whose ancestry (in class A and B stars) is a random variable. 
The underlying physics of this random process is ultimately controlled by the physics of 
r-process nucleosynthesis and the mixing of stellar ejecta. However, even in the absence 
of this detailed physics, simple assumptions allow for a stochastic behavior that resembles 
the observations. We will illustrate two such simple assumptions regarding the random 
distribution of "ancestors." 

First, we treat A^^^, the numbers of ancestors of class A, as a discrete random variable. 
That is, we imagine that the number ancestors from each class take integer values. We choose 
Na from a binomial distribution with mean Na = /o-^Vsn; it then follows that — Ns^ — Na- 
For each star, the fraction of class A (high r) progenitors is — ^a/^SN, and 

R. = fa,.Ra + fb,.Rb ^Rb+ {Ra ' Rb) K/ Ng^ (5) 

The dispersion in Na naturally induces a dispersion in i?*: (7(it!*) = {Ra — Rb) cr{Na)/NsTsi 
Since a{Na) ~ -s/A^, this gives 

a{R,) ~ {Ra - Rb) VK/Ns^ c:^R/^a = {fa Fe/Fe^^^)-'/^ R (6) 

where we have assumed successively that — 1 and Rh ^ R. Thus we see that the fractional 
error in R is just given by the counting statistics of Na'. o'{R^)/R^ ~ l/\/W^. We now can 
put each star on the (r/Fe,[Fe/H]) plane, and the model is complete. 

This model thus has three free parameters which can be fixed by observations: R, Ra, 
and Rh] there is one parameter, (Fe/H)i„in which is (for now) determined by theory. We can 
rewrite eq. (6) explicitly in terms of these parameters: 

<7ii = (/aFe/Fe^i„)-V2^ (7) 

and thus 

InlO 



''As long as both classes A and B exist; below we will consider the effect when there is a time delay before 
the onset of one class. The constancy of r/Fe can be determined observationally. Current data samples the 
low metallicities too sparsely to make a strong statement (Burris et al. (2000)), but with more observations 
one should be able to test for a metallicity dependence of r/Fe. 
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The model thus predicts that the intrinsic scatter in r/Fe falls as (Fe/H)"-*^/^. As we will 
see, this is indeed consistent with the halo star data. Of course, the data also contain 
observational errors, both in [r/Fe] and [Fe/H]. We include these in the model, using 
cr([r/Fe])obs = 0.10 dex, and (T([Fe/H])obs = 0.05 dex, with the supposition that the errors 
are uncorrelated with the intrinsic dispersion. 

Results for a Monte Carlo realization of a population drawn from a binomial distribution 
appear in Figure 1. The predicted envelope is shown, and provides a reasonable fit for the 
observed data. Note that the envelope includes both the observational errors as well as the 
intrinsic scatter given by (cq. 7). In the binomial model just described, there are quantization 
effects due to the integral values of Ni. These are particularly noticeable when Na < 1, which 
occurs when A^sn < l//o — 25, or [Fe/H] < —2.6. In this regime, most points will cluster at 

with only a few at higher values. Furthermore, R — RNa/Nsj^ is also quantized. The 
lower envelope of the nonzero values is given by the case when A^a = 1, so that i?env = -R/A^sn 
and [r/Fe]env = — Fe/H + log^; this is clearly seen in Figure 1. 

One need not require that a star's heavy element ancestry be quantized in this way. For 
example, one can plausibly imagine that the mixing of successive generations of supernova 
ejecta is not an "all or nothing" prospect, but rather that a given parcel of ISM gas and 
dust can be enriched to different degrees by the ancestors it had. This is indeed very likely 
to be the case, which would mean that the quantization effects of the binomial model are 
spurious (and thus the binomial results are only to be trusted in the regime where A^ > 1, 
or [Fe/H] > -2.6. 

Thus, it is of interest to consider making the number of ancestors of each class a con- 
tinuous random variable. This can be done in a way that naturally generalizes from the 
binomial distribution; in that case, the total number of supernova progenitors N^^ is fixed, 
and the fracfAon f = Na/N^-M of progenitors of class A is a random variable which has a 
mean (/) = but can take only discrete rational values due to the integral nature of A^ — a. 
We will continue to take A^sn to be fixed, but we will now assume that the fraction / of 
progenitors of class A is now a continuous random variable, with mean (/) = /„. We thus 
want to draw / from a continuous distribution with values in the interval [0, 1] and with a 
fixed mean. These requirements are met by the beta distribution, which has a distribution 
function 

for / e [0,1] and where B{a,b) is the beta function, which can be expressed in terms 
of the gamma function B(a,6) = r(a)r(6)/r(a + b). We fix the parameters a and b of 
the distribution by simultaneously fitting the mean (/) = a/{a + b) to the observed value 
/a, and forcing the variance (7^(/) = ab{a + b)~'^{a + 6 + 1)~^ to the same value as the 
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binomial case, (J^(/) = /„(! — /a)/-/VsN- These conditions are satisfied ii a — /a(-A^SN — 1) 
and6=(l-/„)(7VsN-l). 

Our procedure is to draw / according to the beta distribution ^(/). We then find 
Na — /A^sN and Nf, — {1 — /)A^sn = -^sn — -^o- Since / is a continuous variable, so are Na 
and Nb. Our choices of the parameters a and b guarantee that (R) — R and that cr{R) is 
identical to the value found for the binomial case, so that eqs. (6) and (7) still hold. Results 
for a Monte Carlo simulation of a stellar population appear in Figure 2. We see that the 
general trend is quite similar to that for the binomial case, but that the continuous nature 
of the parent distribution guarantees that quantization effects are now absent. We see that 
at the very lowest metallicities [Fe/H] < —3.5, the points cluster at r/Fe^^^, and have a 
mean value below that at higher metallicity. This feature has as similar origin to that of its 
counterpart in Figure 1. At the lowest metaUicities, the small number statistics in the SN 
parents is a dominant effect, and the r-rich events are too rare to be expected in a sample 
of the size presented here. 

A comparison with the observed data in Figure 2 shows that the beta distribution model 
gives a good fit to the available data. We thus conclude that the observed scatter in r/Fe 
may be understood by our simple scenario in which halo stars stochastically sample the 
yields of two different populations of heavy element producing events. In this interpretation, 
the degree of scatter increases with decreasing metallicity due to counting statistics, so that 
the lowest-metallicity events record the nucleosynthesis of a few (iVsN = Fe/Femin) events. 
Furthermore, since a{R)/R ~ ^/J^/Nsn, a good measure of the scatter in r/Fe allows one 
to find iVsN and thus one can fix the parameter (Fe/H)min- In this way, one can empirically 
measure the single-event iron yield, and also determine the number of nucleosynthesis events 
recorded at a particular metallicity. The fact that the observed distribution is reasonably 
enclosed by the theory curves which use [Fe/H]niin = —4 (e.g., Audouze & Silk (1996)) 
indicates that this value is a reasonable first approximation. 

2.1.2. Introducing Stellar Lifetimes 

Having achieved a reasonable fit to the data, we are emboldened to consider refinements 
to this model. Specifically, we turn to the question of timescales. Thus far, we have implicitly 
assumed that there is always the same relative probability of r-rich and r-poor nucleosynthesis 
events. However, this assumption is likely to break down at very early times. One expects 
that the to classes A and B of nucleosynthesis events stem from different physics. For 
example, the r-rich and -poor classes could correspond to progenitor stars with different 
masses and thus different lifetimes, or r-process production could arise from binary neutron 



-10- 



star mergers, and thus require some time delay for inspiral. In either case, one expects a 
time lag between the first events of one class versus the first events of the other. During 
this initial period, the r/Fe ratio of the ISM and of any new stars will be only that of the 
allowed class of nucleosynthesis events. There will thus be no dispersion in r/Fe or rather, 
the dispersion will be due to the smaller intrinsic dispersion within the allowed class. Once 
events of the other class can occur, r/Fe dispersion will set in, and the r/Fe will scatter 
around its mean value. 

Wc can crudely simulate this time delay between the onset of the two classes as follows. 
We note that the iron abundance increases with time. Thus, a time delay can be encoded 
as a delay in Fc/H. We choose a value (Fe/II)cut as the cutoff that marks the onset of the 
lagged class of nucleosynthesis events. To allow for dispersion in the birth times of the first 
stars, we introduce some randomness in this cutoff by making the cutoff number A^cut of 
events a poisson random variable, with mean TVcut = Fecut/Femin- We then require that if 
-^SN < -^cut) only events of the allowed class i can occur, and thus R — Ri. If A^sn > -^cut) 
then we follow the standard procedure described above. 

Results appear in Figure 3, using Fecut = 10 Femin, which roughly corresponds to 
timescales riong = 10 Tghort, where the longer lifetime is that of the species cut off. As 
expected, we see that the very low Fe/H, the dispersion vanishes and the points lie at the 
level of the allowed class. While this is clearly an oversimplified description, it nevertheless 
gives a qualitative sense of the kind of behavior one expects if there is indeed a significant 
disparity between the timescales for the production of r-rich and -poor events. A discrepancy 
in timescales is required if the r-process production "switch" is related the progenitor star's 
mass. As we suspect the mass to be the controlling parameter, our prediction is that one of 
the panels in Figure 3 represents the trend that will be observed in ultra-metal-poor stars. 

Thus, a discovery of the sort of behavior seen in Figure 3 would indicate which class has 
the sorter production timescale, and thus provide a key clue as to the astrophysical origin 
of the dominant r-process site. Furthermore, a determination of Fecut/Femin would help to 
quantify (in a model-dependent way) the time delay itself. If class A stars have the longer 
timescale, Ta > r^, then physically this means that (1) the r-process producers evolve more 
slowly and thus (2) these stars should show up later, so that there should be a limiting 
low metallicity below which all stars have i?* = Rh. The appearance of the first stars with 
R^ > Rfj marks the time of the first r-process stars and thus t >Ta- 

We can be more specific by using the fact that r-process stars must be a small fraction 
of all supernovae. This means that the progenitors must represent a limited range in the 
massive star IMF: either 8 — IOM0 or > 25Mq. If is long, this would select the 8 — IOM0 
range, and thus mean that Ta — r(lOMQ) ~ 3 x 10^ yr. Now consider the opposite case. If 
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A stars have the shorter timescale, the Ta < Tb, and the r-process stars should appear first. 
Consequently, there should be a limiting low mctallicity below which all stars have — Ra, 
and the appearance of the first stars with i?^ < Ra marks the time of the first non- r-process 
stars and thus t > Tb- The short implies the high mass (m > 25Mq )progenitor mass 
range, and thus Tb > t{25Mq) ~ 10^ yr. 

While massive stars seem required to explain the r-proccss simply due to their evo- 
lutionary timescale, it is possible that the r-process is not made in the core collapse and 
explosion, but rather through the merger of supernova remnants, via neutron star-neutron 
star mergers (e.g., Eichler, Livio, Piran, & Schramm (1989); Rosswog, Davics, Thielemann, 
& Piran (2000)). As pointed out recently by Qian (2000), an important constraint on this 
scenario comes from magnitude of the halo star r/Fe scatter, which appears to be too small 
(!) to reconcile with the low merger rate and thus high yields of the coalescence events. 
An additional issue for neutron star coalescence is the timescale of the merger, which is an 
extra delay in the appearance of the r-process compared to the other supernova products 
(such as iron). This timescale is a strong function of the orbital semimajor axis and par- 
ticularly eccentricity (Peters (1964)). As shown in population synthesis calculations (e.g.. 
Fryer, Woosley, & Hartmann (1999)), this leads to a wide distribution of inspiral timescales, 
spanning a few Myr to ~ 1 Gyr. Thus, if the r-process does come from binary neutron star 
mergers, then the first events had to have the correct orbital parameters to allow a very 
rapid evolution. 

Finally, even if observations establish the lack of a signature of the kind seen in Figure 
3 would provide useful information about the r-process as well. We would learn that there 
is no significant difference in the timescales for production of the two classes of events. This 
would rule out any models which require such a delay, and would challenge the presumption 
that the r-process yields are a function of the mass of the progenitor stars. Alternatively, it 
may suggest a very different scenario entirely, as we now shall see. 

2.2. r-process-dominated Models 

In the previous section, we assumed that r-process nucleosynthesis events are rare, while 
iron production events are ubiquitous (for massive stars). We now turn to the reverse case, 
which has been developed in some detail by Wasserburg & Qian ((2000); hereafter, WQ2000). 
The essential ingredients are the following, in chronological order: (1) "Pop III." A very rapid 
burst of star formation produces an uneven Galactic fioor of [Fe/H] ~ —4 to —3, but no 
r-process. The inhomogeneity of this population will lead to r/Fe scatter in this metallicity 
range. (2) "Pure r-process." The first of the Pop II supernovae produce the (heavy) r-process 



- 12 - 



without any Fe. These are WQ2000's high-frequency H events. (3) "Fe added." After ~ 10 
i7hfetimes, the first iron-producing supernovae go off. These make the hght r-process (i.e., 
sub-barium, which we do not consider here), but none of the heavy r-process. They occur 
at lower frequency and thus are dubbed L events. 

We model the Wasserburg & Qian scenario as follows; our approach and some details 
are very similar to the model discussed in the previous section. WQ2000 note that the 
Pop III events need not occur in a homogeneous or well-mixed way, which will lead to 
global [Fe/H] variations between —4 and about —3. Thus, we assume that a single Pop 
111 parent would have ejecta which mix with it surroundings, leading to an enrichments of 
[Fe/H] = [Fe/H]mm,iii = —4; this is the Pop III iron yield. This sets the Fe floor. These 
events do not make any r-process, which raises the possibility of an observable abundance 
record of this stage remains, as we will discuss shortly. It is possible, indeed it is likely, that 
the first stars of the next stage of Galactic evolution form prior to the onset of first H events 
which make the first r-process elements. Any sTich stars will have only Pop III abundances 
and thus will be void of any r-process. The number of such r-free stars depends on the 
details of star formation immediately after Pop III, and thus arc difficult to estimate. We 
will thus make the simple assumption that the first H events are not significantly delayed 
after the formation of the first halo stars, so that there is not a significant population of Pop 
III, r-free stars. 

In the next stage, the Pop II stars commence. The purc-r //events occur first, as they 
have lifetimes 10 times shorter than those of the Fe-only L events. Each //event is taken to 
produce the same r-process yield, and thus the total r-production is just the number Nh of 
parents times the yield of a single event. Once the L events commence, they produce only 
Fe (and the hght r-process, which we don't consider). WQ2000 put the Fe yield of a single 
event to be about [Fe/H]^ = —2.5, i.e., 30 times more than that of the Pop III events. We 
assume all L events have this yield. 

We thus have two epochs to model. The first is the period after the Pop III events but 
before the first L event, and the second is the rest of the Pop II phase, up to [Fe/H] = —1. We 
begin modeling each halo star's ancestry by determining which of these epochs it samples. 
To do this, we choose a metallicity in the range —4 < [Fc/H] < — 1.^ We note that whenever 
[Fe/H] < [Fe/H]^ = —2.5, then the star must have arisen in the first, //-only epoch, and have 
no L event progenitors. In this case, the number of Pop HI progenitors is then fixed to be 



^Again, we are free to sample any metallicity in this range, as we are not modeling the massive stars but 
only the present-day halo stars which record the nucleosynthetic past. 
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Nui = Fe/Fe^in,iii.^ 

We then determine the number of i? events in our halo star's hneage, by choosing Nh. 
As in the previous section, we note that the specific distribution of is determined by the 
detailed physics of star formation and gas dynamics in the earliest epoch in our Galaxy. The 
distribution is constrained by the H event-to- L event ratio of 10, which demands that on 
average, there are 10 H events before the first L event occurs. Independent of the specific 
distribution, a crucial issue is the expected correlation between and Nm (and thus Fe). 
Since these events occur in two different phases by two different population, we will make 
the simple assumption that the two are uncorrelated. We then have i?* = rhNh/Fe, with Fe 
already fixed and with Vh the r-process yield of a single H event (determined below). Note 
that since is uncorrelated with Fe, then the mean trend (R) — {N^jrh/Fe which scales 
as Fe"^. That is, the mean r/Fe trend for —4 < [Fe/H] < —2.5 should not be constant, but 
should decrease with logarithmic slope -1. Scatter about the r/Fe trend is created by the 
variance in Nh. = a{Nh)rh/Fe = a{Nh)/Nh R. For illustration, we will choose from 
a uniform distribution, in (0, 10), which gives {R) = 5rh/Fe, and gr/R = 1/^3. 

We now turn to the case of stars for which we have picked a higher Fe/H > Fe/H^. 
For these stars, some contribution by iron-producing L events is necessary. As before we 
generate a number iVsN of supernova parents for each halo star Thus, each star generated 
has Nh H parents and A^^ L parents, but the star's iron abundance constrains our choice of 
A^ and thus A^sn- Namely, the star has a metallicity 

(Fe/H)^ - Ari„(Fe/H)„i„,„i + Ar,(Fe/H), = {Nm,eS + N,) (Fe/H), (10) 

with Aiii.cff = (Fejnm,iii/Fe,) Aiii. We choose A^m from a distribution between and Fe^/Feminjn ~ 
30; for illustration, we will use the uniform distribution Once given Ahi, then A, is fixed by 
eq. (10). 

Turning to N^, we must choose this in a way that is consistent with our choice of A,. 
Note that since the L events are 10 times rarer, the variance in A'^ is much larger than that 
in Nh, and thus we should expect a large span in the allowed Nh choice at fixed A,. We 
proceed by first estimating the total number of supernova events to be iVsN = A,//,, where 
f I = 1/11 — 1 — fh IS the fraction of supernova of type L. Using a beta distribution, we 
choose /,,^, with a = ft{N^^ - 1) and h = ft{Ns^ - 1). We then have = 1 - and 
can determine A^sn,* = A^^//^,*- From these, we find A^^ = fh,*NsN,* = fh,*/fe N^, and the 



^This may not be an integer, but we will interpret fractional progenitor numbers as indicative of partial 
mixing, as done in the previous section. 
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star's r/Fe ratio, 

where i?yid = Th/Fei is the ratio of the (number) yields from the two event types. Note that 
at high metaUicity, we have N^, 3> iVin,eff and thus 

R* — -Ryid — R (12) 

But since Wasserburg & Qian specify that H events occur 10 times more frequently than 
L events, we have Nh/Ni = 10, and from eq. (12) we can infer the proper Ry\d = Vh/Fee 
parameter to give the correct R. Namely, we have i?yid = R/10 = 0.2, and r/j = FeiRyu = 
6.3 X 10-1 

We now turn to the r/Fe dispersion at Fe/H > Fe/H^. Note that simple analytic formu- 
lae are not easily derived for the r/Fe scatter in the WQ2000 scenario, owing to the presence 
of random variables in both the numerator and denominator of eq. (11). Nevertheless, one 
can roughly estimate that a{R)/R ~ l/y/W^. For Fe > Fe^, this gives an intrinsic dispersion 
of a{R)/R ~ y/Fe/Fee, the same Fe^^^ scahng as that in the Fe-dominated case. 

Results from a Monte Carlo realization appear in Figure 4. the dashed curve shows 
the 2(7 envelope which includes observational errors and the intrinsic scatter predicted by 
a{R)/R = ^/Fe/Fe^ and plotted over its range of vahdity. We see that the scatter at 
[Fe/H] > —2.5 is similar to that in the Fe-dominated case, and in reasonable agreement 
with the present data. However, the low metaUicity convergence behavior seen in Figure 3 
does not appear in the Wasserburg & Qian scenario. Rather, at [Fe/H] < —2.5 we see the 
R ~ 1/Fe trend discussed above. We again conclude that observations of the r/Fe scatter 
in this regime can thus provide important information about r-process astrophysics. 



3. Observational Tests of r-Process Astrophysics 

The models we have presented for r/Fe scatter all show similar trends from Fe/H = — 1 
down to about —2.5. However, below this metaUicity, the differences become increasingly 
stark, as seen in the figures. The predicted behaviors are for points to cluster above, below, 
or around the mean value [r/Fe]avg- 

• If the points all cluster around [r/Fe] ^ [r/Fe]avg, this points to a very early population 
which, regardless of the small number of ancestors, maintained high r/Fe. This would 
confirm that the earliest stars produced both iron and r-elements, and thus that the 
sources of the r-process are short-lived and thus high-mass stars. 
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• If, conversely, the points all cluster around [r/Fe] <S [r/Fe]avg (or which have only 
upper limits for r-process elements), this indicates that the earliest stars produced iron 
but failed to produce significant r. This is thus evidence for a delay until the first 
r-producing stars, and thus points to longer-lived progenitors (lower mass supernovae 
or NS-NS binaries). 

• Finally, the points could continue to cluster around [r/Fe]avg- This would indicate 

that the earliest stars could collectively produce both the r-process or iron, but the 
persistent scatter would demand that individual stars could only produce one of the 
two. This points to a picture of the kind suggested by Wasserburg & Qian, wherein 
an early (and non-uniform) Pop 111 production of iron is quickly followed by rapid 
production of iron-free r-process ejecta. 

Thus we see that ultra-metal-poor stars can provide "smoking gun" evidence for r- 
process astrophysics. 

4. Conclusions 

We have presented an analytical approach to understanding the scatter in r-process- 
to-iron ratios in metal-poor halo stars. The observed r/Fe scatter demands that the bulk 
of the r-process and of iron cannot be produced in the same stars. We have used a simple 
stochastic description of the different populations of nucleosynthesis events which contribute 
to r-process and iron abundances in each halo star. The random star-to-star variations in 
nucleosynthetic "ancestry" lead to scatter in r/Fe. The models we present are all successful 
in reproducing the scatter in the available data, which go down to about [Fe/H] = —3. We 
see that the abundance scatter can be understood simply in terms of the variance due to 
counting statistics of the ancestry of each star. We present simple and convenient analytic 
expressions for the scatter in this regime. 

We model two basic scenarios for r-process production; one in which the r-process is 
rare, occurring in about ~ 4% of all massive stars. The other scenario, that of WQ2000, takes 
the reverse approach in which r-process production is ubiquitous and iron production rare 
in massive stars. ^ We find that these scenarios give very similar r/Fe scatter at [Fe/H] > 
—2.5, but show divergent behavior at lower metallicity. Observations in this regime can 
discriminate among these scenarios and give important clues regarding the lifetimes and 



''In the final stages of completion of this paper, we became aware of the work of Qian (2001), which also 
considers r/Fe scatter in the WQ2000 model. His methods and conclusions are very similar to ours. 
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thus the masses of the astrophysical r-process production sites. We therefore strongly urge 
an observational effort to measure r-process abundances (and upper limits!) in the ultra- 
metal-poor halo stars. 

We thank Roberto Gallino, John Lattanzio, Chris Sneden, and Ron Webbink for useful 
discussions. This research has received support from NSF grant AST-9986974 to J.J.C. and 
from DOE contract B341495 to J.W.T. 
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Fig. 1. — Distribution of r/Fe abundances for halo stars. Upper panel: model calculation. 
The points represent stars whose ancestry has been drawn according to a binomial distribu- 
tion. In addition to the intrinsic dispersion calculated in the model, we have included random 
observational errors with zero mean and standard deviation cr[/j] = 0.15 and crjpe/H] = 0.10. 
Lower panel: observed data for [Eu/Fe]; references are as shown. For CS 31082-001 (Cayrel 
et al. (2001)) Eu abundances are as yet unavailable, so the average of [Os/Fe] and [Ir/Fe] is 
used. The dashed curves in both panels are the analytic relation for the 2 — a envelope pre- 
dicted from this model (eq. 7). The curves are drawn to account for the logarithmic nature 
of the plot, ie., the curves are located at the 2a limits, namely log^ ± 2\ogl0a{R)/R. 
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Fig. 2. — As in Figure 1, for the beta distribution model. Due to the continuous nature of 
the parent distribution, the quantization effects seen in the binomial case are now absent. 
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Fig. 3. — Predicted r/Fe distributions, with timescale information encoded through metal- 
hcity cutoffs. Upper panel: a delay has been imposed on class B (r-poor) events, so that 
class A events dominate at early times. Lower panel: a delay has been imposed on class A 
events. 
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Fig. 4. — As in Figure 1, for the Wasserburg & Qian (2000) model. 



